Coexistence of wrinkles and blisters in supported graphene

Blisters induced by gas trapped in the interstitial space between supported graphene and the substrate are commonly observed. These blisters are often quasi-spherical with a circular rim, but polygonal blisters are also common and coexist with wrinkles emanating from their vertices. Here, we show that these different blister morphologies can be understood mechanically in terms of free energy minimization of the supported graphene sheet for a given mass of trapped gas and for a given lateral strain. Using a nonlinear continuum model for supported graphene closely reproducing experimental images of blisters, we build a morphological diagram as a function of strain and trapped mass. We show that the transition from quasi-spherical to polygonal of blisters as compressive strain is increased is a process of stretching energy relaxation and focusing, as many other crumpling events in thin sheets. Furthermore, to characterize this transition, we theoretically examine the onset of nucleation of short wrinkles in the periphery of a quasi-spherical blister. Our results are experimentally testable and provide a framework to control complex out-of-plane motifs in supported graphene combining blisters and wrinkles for strain engineering of graphene.Peer ReviewedPostprint (published version

Hydraulic fracturing in cells and tissues: fracking meets cell biology

The animal body is fundamentally made of water. A small fraction of this water is freely flowing in blood and lymph, but most of it is trapped in hydrogels such as the extracellular matrix (ECM), the cytoskeleton, and chromatin. Besides providing a medium for biological molecules to diffuse, water trapped in hydrogels plays a fundamental mechanical role. This role is well captured by the theory of poroelasticity, which explains how any deformation applied to a hydrogel causes pressure gradients and water flows, much like compressing a sponge squeezes water out of it. Here we review recent evidence that poroelastic pressures and flows can fracture essential biological barriers such as the nuclear envelope, the cellular cortex, and epithelial layers. This type of fracture is known in engineering literature as hydraulic fracturing or "fracking"Peer ReviewedPostprint (author's final draft

An atomistic-based foliation model for multilayer graphene materials and nanotubes

We present a three-dimensional continuum model for layered crystalline materials made out of weakly interacting two-dimensional crystalline sheets. We specialize the model to multilayer graphene materials, including multi-walled carbon nanotubes (MWCNTs). We view the material as a foliation, partitioning of space into a continuous stack of leaves, thus loosing track of the location of the individual graphene layers. The constitutive model for the bulk is derived from the atomistic interactions by appropriate kinematic assumptions, adapted to the foliation structure and mechanics. In particular, the elastic energy along the leaves of the foliation results from the bonded interactions, while the interaction energy between the walls, resulting from van der Waals forces, is parametrized with a stretch transversal to the foliation. The resulting theory is distinct from conventional anisotropic models, and can be readily discretized with finite elements. The discretization is not tied to the individual walls and allows us to coarse-grain the system in all directions. Furthermore, the evaluation of the non-bonded interactions becomes local. We test the accuracy of the foliation model against a previously proposed atomistic-based continuum model that explicitly describes each and every wall. We find that the new model is very efficient and accurate. Furthermore, it allows us to rationalize the rippling deformation modes characteristic of thick MWCNTs, highlighting the role of the van der Waals forces and the sliding between the walls. By exercising the model with very large systems of hollow MWCNTs and suspended multilayer graphene, containing up to 109 atoms, we find new complex post-buckling deformation patterns

Approximation of tensor fields on surfaces of arbitrary topology based on local Monge parametrizations

We introduce a new method, the Local Monge Parametrizations (LMP) method, to approximate tensor fields on general surfaces given by a collection of local parametrizations, e.g.~as in finite element or NURBS surface representations. Our goal is to use this method to solve numerically tensor-valued partial differential equations (PDE) on surfaces. Previous methods use scalar potentials to numerically describe vector fields on surfaces, at the expense of requiring higher-order derivatives of the approximated fields and limited to simply connected surfaces, or represent tangential tensor fields as tensor fields in 3D subjected to constraints, thus increasing the essential number of degrees of freedom. In contrast, the LMP method uses an optimal number of degrees of freedom to represent a tensor, is general with regards to the topology of the surface, and does not increase the order of the PDEs governing the tensor fields. The main idea is to construct maps between the element parametrizations and a local Monge parametrization around each node. We test the LMP method by approximating in a least-squares sense different vector and tensor fields on simply connected and genus-1 surfaces. Furthermore, we apply the LMP method to two physical models on surfaces, involving a tension-driven flow (vector-valued PDE) and nematic ordering (tensor-valued PDE). The LMP method thus solves the long-standing problem of the interpolation of tensors on general surfaces with an optimal number of degrees of freedom.Comment: 16 pages, 6 figure

Blending isogeometric analysis and local maximum entropy meshfree approximants

We present a method to blend local maximum entropy (LME) meshfree approximants and isogeometric analysis. The coupling strategy exploits the optimization program behind LME approximation, treats isogeometric and LME basis functions on an equal footing in the reproducibility constraints, but views the former as data in the constrained minimization. The resulting scheme exploits the best features and overcomes the main drawbacks of each of these approximants. Indeed, it preserves the high fidelity boundary representation (exact CAD geometry) of isogeometric analysis, out of reach for bare meshfree methods, and easily handles volume discretization and unstructured grids with possibly local refinement, while maintaining the smoothness and non-negativity of the basis functions. We implement the method with B-Splines in two dimensions, but the procedure carries over to higher spatial dimensions or to other non-negative approximants such as NURBS or subdivision schemes. The performance of the method is illustrated with the heat equation, and linear and nonlinear elasticity. The ability of the proposed method to impose directly essential boundary conditions in non-convex domains, and to deal with unstructured grids and local refinement in domains of complex geometry and topology is highlighted by the numerical examples

Rippling and a phase-transforming mesoscopic model for multiwalled carbon nanotubes

We propose to model thick multiwalled carbon nanotubes as beams with non-convex curvature energy. Such models develop stressed phase mixtures composed of smoothly bent sections and rippled sections. This model is motivated by experimental observations and large-scale atomistic-based simulations. The model is analyzed, validated against large-scale simulations, and exercised in examples of interest. It is shown that modelling MWCNTs as linear elastic beams can result in poor approximations that overestimate the elastic restoring force considerably, particularly for thick tubes. In contrast, the proposed model produces very accurate predictions both of the restoring force and of the phase pattern. The size effect in the bending response of MWCNTs is also discussed